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ABSTRACT 

This paper explores the potential for statistical epoch of reionization (EOR) measurements using 
wide field radio observations. New developments in low frequency radio instrumentation and signal 
processing allow very sensitive EOR measurements, and the analysis techniques enabled by these 
advances offer natural ways of separating the EOR signal from the residual foreground emission. 
This paper introduces the enabling technologies and proposes an analysis technique designed to make 
optimal use of the capabilities of next generation low frequency radio arrays. The observations we 
propose can directly observe the power spectrum of the EOR using relatively short observations, and 
are significantly more sensitive than other techniques which have been discussed in the literature. For 
example, in the absence of foreground contamination the measurements we propose would produce 
five 3-sigma power spectrum points in 100 hours of observation with only 4 MHz bandwidth with 
LOFAR for simple models of the high redshift 21cm emission. The challenge of residual foreground 
removal may be addressed by the symmetries in the three-dimensional (two spatial frequencies and 
radiofrequency) radio interferometric data. These symmetries naturally separate the EOR signal from 
most classes of residual un-subtracted foreground contamination, including all foreground continuum 
sources and radio line emission from the Milky Way. 
Subject headings: 



1. INTRODUCTION 

The possibility of observing 21 cm emission from 
the "cosmic dark a ges" was first recognized by 
Sunvaev & Zeldovich (1972) and later developed b y 
Scott fc Rees ( 1 99(rL Madau. Meiksin fc R.eesI l|1997D . 
and iTozzi et all l|200f]l) to show the nature of the ex- 
pected signal and demonstrate that future low fre- 
quency radio telescopes could measure this emis- 
sion. This has inspired a significant amount of work, 
with important contr i butions by Ilicv ct ah ( 2002 ) , 
Pi Matteo et alJ i2002tl.lCarini. Gnedin fc Ow^ |2002|i . 
Fiir1a,Tietto fc T.oeH f200ir and others. Just recently, 
Barkana fc Loebl (|2004) showed that the predicted global 
step in reshifted 21 cm emission was an artifact of 
the small box size used in the simulation procedures. 
Combined with the high Th omson scattering cross 
section observed by WMAP l|Kogut et all l2003() and 
other arguments, these results suggest that reioniza- 
tion is much more complicated than originally conceived 
(flui & Haiman 2003). 

These theoretical developments coincide with the de- 
velopment of several advanced low frequency radio ar- 
rays which will come online by the end of the decade. 
Recent advances in digital electronics are enabling the 
construction of radio arrays which directly sample the 
incoming wavefront with thousands of simple antennas 
and perform all signal processing in the digital domain. 
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The digital processing simplifies the mechanical systems, 
allowing very large collecting areas at modest cost and 
enabling entirely new observing modes. 

One of the enabling digital processing technologies is 
the Wide Field Correlator ( WFC) developed at M IT and 
the Haystack Observatory ('Mo rales etaP 12004V The 
digitized signal from up to several thousand individual 
antennas can be fully correlated by the WFC, producing 
more than 1 billion individual visibility measurements. 
These visibilities allow the full field of view of the anten- 
nas to be imaged, and greatly increases the sensitivity 
of epoch of reionization (EOR) observations by allowing 
both parallel observation of a large number of sky loca- 
tions and by enabling new statistical analysis techniques 
(see Section 

In this paper we show how wide field radio observa- 
tions combined with statistical analysis techniques simi- 
lar to those used in CMB experiments may be applied to 
studies of the EOR. An analysis usin g statistical CMB 
techni ques independently developed bv lZaldarriaga et all 
l)2004t) has been recently published, where they show the 
usefulness of the EOR power spectrum and that fore- 
ground contamination can be easily removed. Our ap- 
proach is similar but emphasizes the importance of the 
full three dimensional statistics inherent in the EOR sig- 
nal. By making full use of all the available informa- 
tion, our technique is both more sensitive to the EOR 
signal and provides additional handles for removing the 
foreground contaminatio n. This paper complements the 
iZaldarriaga et ahl l)2004j) paper and details the experi- 
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mental and analysis advances which will enable future 
radio arrays to make robust EOR measurements using 
reasonable integration times. 

This discourse is divided into three major themes. In 
Section 121 the basic design and wide field capabilities of 
future low frequency radio arrays are presented, followed 
by a description of the analysis technique enabled by the 
wide field of view in section |3| Sections 0] and |5l then 
present sensitivity calculations for wide field radio EOR 
observations and the constraints they can provide. Ad- 
ditional calculations and issues important to the EOR 
field are detailed in the appendices. 

2. WIDE FIELD RADIO OBSERVATIONS 

In contemporary and planned radio telescopes, the dig- 
itization happens as early as possible in the signal collec- 
tion process and signal processing is performed in high 
speed digital electronics. This shift to early digitization 
dramatically affects the "look" of the array, with thou- 
sands of small elements replacing a few large parabolic 
antennas. This basic change in the design philosophy en- 
ables wide field observations but also places much greater 
demands on the backend digital processing. This section 
introduces the basic antenna design and the massively 
interconnected correlators that can scale to thousands 
of antennas. While none of the proposed arrays have 
been constructed, they all share the same basic design 
elements which are used as the basis for the following 
discussion. 

The individual receptor elements consist of simple 
dipoles. Dipoles may then be grouped with an analog 
beamformer into antenna "tiles." The more dipoles per 
tile, the larger the collecting area and smaller the field 
of view. Depending on the the science goals, the num- 
ber of dipoles per antenna is then chosen to optimize 
science return and cost. The 80 - 240 MHz prototype 
tile developed by the MIT Haystack Observatory con- 
sist of 16 crossed dipoles, maximizing the collecting area 
at frequencies critical for EOR studies at the expense of 
limiting the field of view to 20° (see Figure ^ . After 
each antenna the signal is immediately digitized and all 
subsequent processing is done in the digital domain. 

Wide field observations which image the full antenna 
field of vie w are enabled by the Wide Field Correlator 
technology (! Morales fc Cappalloll200^ . Since the num- 
ber of visibilities grows as the square of the number 
of antennas, arrays with thousands of antennas require 
millions of baselines, necessitating a massively intercon- 
nected correlator technology. The WFC technology en- 
ables large arrays to be fully correlated at low cost by 
developing scalable techniques for handling millions of 
pairwise connections with a minimum of point-to-point 
connections. The WFC technology also addresses a num- 
ber of additional technological issues to allow efficient 
use of FPGA chips including time multiplexing to match 
bandwidth with chip speed, and data buffering and re- 
ordering to minimize on chip memory requirements. This 
architecture is very scalable, and the performance of fu- 
ture versions of the WFC should increase with Moore's 
law providing even greater performance at very reason- 
able cost. 

The first generation wide field correlator (WFCl) per- 
forms the pairwise correlations for 3072 antennas at 4 
MHz bandwidth, producing 1.2 billion 16 kHz visibilities 



fsee lMorales et al.) l)2004j) for a full description). The ca- 
pacity of an FX correlator design can be approximated 
by the trillions of complex multiply and accumulates 
per second (Tcmacs) performed, given by the number of 
baselines (pairwise combinations of antennas) times the 
bandwidth. The WFCl has a performance of 38 Tcmacs 
( (30722/2) basefines x 4 MHz x 2 pol = 38 Tcmacs), 
and is expected use ~ 2700 FPGA chips on 42 boards 
with a hardware cost of under 1 million dollars US (2006 
price). For massively interconnected correlators such as 
the WFCl, the number of antennas can be traded for 
bandwidth with little additional cost as long as the num- 
ber of Tcmacs remains constant. 

The baseline observing program for wide field EOR 
observations involves a few hundred hours of observa- 
tions in a cold portion of the sky. While all of the fore- 
ground sources will be subtracted from the data, residual 
foreground contamination can be reduced by carefully 
choosing the observation direction. The observation lo- 
cation will be away from the galactic plane to reduce the 
strength of the galactic emission, and will be chosen to 
contain a minimum of bright sources in an area where the 
galactic emission is particularly smooth. The observing 
times are selected when the chosen location is at high 
elevation angle and the ionospheric conditions are par- 
ticularly favorable. In practice this will mean observing 
for a couple of hours on many nearly consecutive nights 
when ionospheric conditions are at their best. With six 
months of observing, a very clean data set can be ob- 
tained with a deep integration of a few hundred hours. 

The combination of very large effective area and wide 
field of view makes next generation low frequency obser- 
vatories ideally suited to EOR measurements. While the 
wide field of view can be thought of in the traditional 
sense as allowing ~ 1 hundred thousand simultaneous 
pointings, in the context of EOR observations it is more 
productive to think in terms of the visibility measure- 
ments. As will be presented in the next section, the EOR 
information is most easily extracted in the Fourier repre- 
sentation, and one can think of the visibilities produced 
by the WFCl as 1.2 billion independent measurements. 
This large number of independent measurements allows 
wide field radio observations to directly sample the sta- 
tistical signature of the EOR. 

3. STATISTICAL ANALYSIS TECHNIQUE 

In this section we introduce our statistical analysis of 
EOR measurements in three parts. In the first part we 
introduce the formalism for calculating the three dimen- 
sional statistics of the EOR signal. This is followed by 
a discussion of the inherent symmetries of the statis- 
tics and how this can be used to separate residual fore- 
ground contamination from the EOR signal in Section 
13.21 and a discussion of the non-Gaussian properties of 
the statistics in Section 13.31 The formalism presented 
in Section |3.1l parallels calculations for the CMB by 
iWhite et al.l l|1999ll extend ed to three dimens i ons, a nd is 
similar to calculations bv IZaldarriaea et all lj2004l) and 
IBharadwai fc Sethil lj200H) but uses Fourier formalism 
to clearly illustrate the three dimensional nature of the 
problem and to highlight the symmetries and statistics 
which are explored in Sections 13.21 and 13.31 The idea 
of using power spectra to measure the EOR signal has 
been around for some time, but the power of the symme- 



Fig. 1. — This photograph shows a prototype for one sixteen element antenna for the 80 - 240 MHz band. 



tries inherent in the three dimensional data set and the 
non-Gaussian statistical properties of the signal have not 
been previously appreciated. 

3.1. Introduction to the Visibility Statistics of the EOR 

In wide field radio observations, the observed frequency 
of the redshifted 21 cm line maps to the line of sight dis- 
tance of the source. While peculiar velocities mean this 
mapping is not exact — a small effect to be studied in fu- 
ture work — an approximate image cube can be formed by 
replacing the observed frequency with the equivalent line 
of sight distance. The resulting image cube is a full three 
dimensional tomograph of neutral hydrogen emission in 
the high redshift universe. 

The emission strength from neutral hydrogen during 
the epoch of reionization depends on the local hydrogen 
density, the fractional ionization, and the spin temper- 
ature of the gas, and gives us a differential specific in- 
tensity A/ at the observed location with respect to the 
cosmic microwave background radiation (see Appendix 
^for a full derivation of A/). Because the 21 cm line is 
narrow and weak and the gas is optically thin, there is no 
significant absorption of background emission by neutral 
hydrogen nearer the observer. This allows us to directly 
form an image cube which represents the three dimen- 
sional emission of the neutral hydrogen. Furthermore, 
a three dimensional Fourier transform can be applied to 
the image cube to determine the spatial structure of the 



EOR signal. This gives us four representations of the 
signal: 

AI{f) Ai{k) 



(1) 



AI{6^,6y,Af) ^ AIiu,v,7j). 
The horizontal pairs in Equation^are related by Fourier 
transforms and the vertical pairs are related by a change 
of variables from the observational quantities sky posi- 
tion Ox,Oy and frequency A/ (and their inverses u,v,ri) 

to the comoving distance r (and inverse k). The upper 
left hand term is the brightness distribution of the 21 
cm emission at the source while the lower left hand term 
is the observed image cube. Similarly the upper right 
hand term gives the spatial structure of the EOR signal 
and the lower right hand term is the spatial structure 
of the three dimensional image. The mapping between 
variables is given by 



A/: 



y dm(z)- 

_ Hof2iE{z) A 



27r ' 

kyPMjz) 

c(l + z)^ 

2TTHof2lEiz) 



(2) 
(3) 
(4) 
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where /21 is the frequency of the 21 cm hn e, Dm is 
the transverse comoving distance (jHogdllOOflj) . E{z) = 

y/flMC^- + z)^ + f^k(l + z)^ + ^A, and z is the reference 
redshift at the middle of the observed data cube used for 
determining A/ and Ar^.^'^ 

For an interferometric array the fundamental observ- 
able is the set of cross-correlations between the antenna 
elements V{u,v, Af), called visibilities in radio astron- 
omy, and is related to the image cube and the Fourier 
representation via Fourier transforms. The Fourier rep- 
resentation of the image cube AI{u,v,r]) is related to 
the measured visibilities via the one dimensional Fourier 
transform along the frequency axis 



(10) 



and the image cube is obtained by performing a two di- 
mensional Fourier transform in the sky coordinates 



Aiie,,ey,Af) 



V{u, V, A/)e'27r(«e.i+t,e„j)^^ 



These relationships are shown graphically in Figure [3 

In order to calculate the actual measured values of A/, 
A/, and V, we need to understand the mapping from the 
true intensity distribution AIt to the measured quanti- 
ties. In particular, the finite field of view and bandwidth 
of the measurement truncates the Fourier transform from 

^ The top line of Equation is a familiar relationship from cos- 
mology and the bottom line is a standard of radio astronomy. Un- 
fortunately the two different communities use different conventions 
for the Fourier transform. In this paper we will use the transfor- 
mation pair 

A/{r) = / /\i{k)e^^-^d?k, 



AI{k) 



(27r)3 



(5) 
(6) 



for transforming functions of position f and k and 

(7) 



AI(u, V, rj) 



Ai(e^,ey, A/)e-'2'^("''-'+"''!'j+'"^^'^'de:, dOy dAf, 

(8) 



for functions of the observation variables. While using two differ- 
ent definitions of the Fourier transform is awkward, it allows us 
to use the standard notation of both cosmology and radio astron- 
omy and avoids the inevitable confusion that would result from 
non-standard notation. It is this difference in transforms which 
explains the relative 2it factors in Equations TUVil Standard units 
are comoving Mpc for r, rad per comoving Mpc for k, radians for 
6, Hz for A/, wavelengths (or equivalently rad~^) for u and v, and 
Hz~^ for r^. If r and k are in proper coordinates instead of comov- 
ing coordinates. Eg nations l2ll4l need to be divided by (I-I-2). Note 
that while r] physically has the units of seconds, it is really acting 
as a stand-in for the spatial structure and is best thought of as an 
inverse distance. 

^ An exact derivation maps the observed frequency / to 



Ho Jf^^^ f2iDM{f')E(f') 



(9) 
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is then an exactly linear map to the line of sight distance and has 
e same scale as 9^ and 9y. This alternate development leads to a 
kernel which must be adde d to the Fourier transform relationship 
developed in Eouation 1101 and Figure |2] For almost all reasonable 
parameters the approximate relationship in Equation|l]is sufficient. 



the true brightness distribution to the measured visibili- 
ties. 

In what follows we neglect the curvature of the sky. For 
the purposes of examining the statistics of radio emis- 
sion on the celestial sphere, this assumptions is valid 
when 2sin{0/2) « 9 l|Bond fc Efstathiou 19871 For our 
proposed observations the largest value of 9 is 10°, and 
we can adopt the two dimensional Fourier relationships 
rather than the full statistical treatment of spherical har- 
monics. The condensed arrays are expected to be copla- 
nar to a fraction of a wavelength at FOR frequencies, 
so the "w-term" from relative antenna elevations may be 
neglected. Additionally, the visibilities from the Wide 
Field Correlator will be collected every 1/2 second and 
the baseline adjusted to allow tracking of a single patch 
of the radio sky. 

The true distribution of the FOR specific intensity is 
given by 

^™ (12) 



A/t 



HI 



C 



(ph)' 



where Fhi gives the fluctuations in the neutral hydrogen 
emission. This can be separated into a spin temperature 
component C and the density terms, where pn and pm 
are the density of hydrogen and neutral hydrogen respec- 
tively. For a baryon density JIb/j^ = 0.02 and a helium 
fraction Y ~ 0.24, C is given by 



Cjy = (2.9 mK) h 
when measured in Jy str~^, and 



_i (-fffl — -^cmb) 



E{z) Ts 



(13) 



Ck ^ (2.9 mK) h-' ^^+f (Ts-Tcmb) ^^^^ 



in brightness temperature units. C captures the spin 
temperature dependence of the neutral hydrogen emis- 
sion and Phi/ (ph) includes the effects of density fluctu- 
ations and ionization fraction, and in general both terms 
are functions of position.'^ A full derivation of C along 
with definition of all terms is given in Appendix 1X1 

The first step in determining the measured quantity 
is to define the window function W{6x,9y, f) given by 
the instrumental field of view and bandwidth. This can 
be Fourier transformed to obtain W{u,v,r]). Since the 
measured brightness AI{9x,9y, f) is given by multiplying 
the true brightness by the spatial window function, the 
measured Fourier transform of the intensity is given by 
convolving the Fourier transform of the true intensity 
A/t(u,w,?7) with W{u,v,r]) 



Ai{ui,vi,rii) = 
AIt{u, V, r])W{ui — u,Vi — v, rji — rf) du dv dr/. (15) 



^ We have followed EOR convention and defined AIt with re- 
spect to the CMB radiation and not the CMB radiation plus the 
average contribution of the neutral hydrogen emission. This choice 
of zero point differs from the convention used in derivations of 
CMB fluctuations, and can add a <5-function at zero wavenumber 
to the power spectrum. This offset is very hard to measure in 
practice, and its contribution is ignored in the remainder of this 
paper. Mathematically the zero point contribution can be removed 
by subtracting C (pm) / (ph) from Equation 1121 to redefine AIt- 
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Fig. 2. — This diagram shows the Fourier transform relationships between the image cube, the measured visibiUties, and the Fourier 
representation. For an interferometer, the fundamental observable is the visibility-frequency cube, which can then be transformed into 
either an image cube to show source locations or the Fourier representation to analyze the spatial structure of the signal. 



Using Equation [T^ this can be rewritten as 

A/(ui,wi,77i) = 

Fhi (u, V, ri)W{ui — u,vi — v, rji — rf) du dv drj (16) 



for a specific realization of Fhi. Because model predic- 
tions of Fm are statistical it is convenient to consider 
the average value. For the next several steps we will first 
sketch the derivation of the statistical relationship be- 
tween the intensity and the neutral hydrogen fluctuations 
in n space, where n is a wave number vector in undefined 
units, but is used to stand in for k or ui + vj + r]k. After 
working through the derivation in the abstract n space 
we will redo the analysis in physical units. Rewriting 
Eauation lltil in the n space we obtain 



A/(n') = I Fm{n)W{n' ~n)d^n 



(17) 



We can use this to calculate the expectation value of two 
visibilities measured at two different locations n' and n" 



Aiin') Ar{n") 




(rt) W{n' - n) d^i 



F^ln"')W*{n" -n"')d^n" 



(18) 



The model expectation for the fluctuations is spatially 
homogeneous (no preferred point in space), so the expec- 
tation for the Fourier transformed value Fm (n) has ran- 
dom phase. The statistical average (Fmin) F^-^{n'")) — 
(|FHi(n"')P) S{ri - n'") d^n'" due to the random phase, 
and this is not changed by multiplying by a deterministic 
function like W. So the expectation value gives us: 

A/(7l') Ar(n")) = J Pmin) W{n'-n) W*{n"^ri) d^n, 

(19) 

where Pm{n) = (|FHi(ri"')p) 5{n — n"') and is the power 
spectrum of neutral hydrogen emission.'* This is the re- 
sult we were looking for which relates the statistical prop- 
erties of the intensity to the power spectrum of emission 
fluctuations, and is equivalent to the standard CMB de- 
velopment extended to three dimensions. For cases of 

Note that sometimes the power spectrum is defined as 
Phi(") = (27r)-3 (|Fhi(*i"')P) - n'") in the CMB literature 
due to a different definition of the Fourier transform relation (see 
Footnote [Tl Equations ISl and ISt . leading to a difference of {2-k)~^ 
in the results. 



constant ionization fraction and spin temperature, the 
analogy is even closer and Phi becomes simply the power 
spectrum of the matter density fluctuations (see Section 
[S)) . This equation also shows a number of interesting fea- 
tures including the correlation length of the intensities. 
The easiest way to use this in an analysis is to compare 
the intensity of locations with themselves to obtain the 
power spectrum of the intensity fluctuations 



A/(n') ) = Pm{ri)\W{ri' -n)\ 



'd^n. 



(20) 



We can now repeat these steps using the observation 
units u,v,ri. Eauation 1191 gives the expectation for the 
intensity at two locations: 



A/(ui,wi,?7i) Ar{u2,V2,r]2] 



Phi(m, V, v) W{ui ~ u,vi — V, rji — ri)x 

W*{u2 — u,V2 — V, 772 — 77) du dv drj. (21) 

For a single location this simplifies to 

~ 1 2 \ 

A/(ui,i;i,77i)| 

Pm{u, V, 77) |iy(ui — u,vi — v,r]i — ri)\^ du dv drj. 

(22) 

For the specific case of observing a single small area of the 
sky, Parseval's relation can be used to convert the con- 
volution in Equation 1221 to the integrals over the Fourier 
space found in Equ ations 1 2 and 13 of Tozzi et al. (2000), 
as corrected bv lllliev et al.l (|2p02). 

The statistical distribution calculated in Eauation 1221 
can be directly sampled using the 1.2 billion visibilities 
produced by the first generation Wide Field Correlator. 
Because of the correlation length introduced by the win- 
dow function, not all of the visibilities measure inde- 
pendent values of A/ (there are ~ 20 million indepen- 
dent values of A/). However, because the visibilities are 
independent measurements their noise is uncorrelated,^ 
and the uncertainty on values of A J which are measured 
by multiple visibilities is correspondingly reduced. The 
large number of well constrained A/ measurements en- 
abled by wide field radio observations allow the statistical 
distribution in Equation|221to be directly sampled. In the 

^ There is some evidence of correlated power on short baselines 
from the VLA. From a theoretical point of view this must be either 
an instrumental effect or a foreground signal such as the clustered 
faint galaxy population. In either case this is an ef fect which must 
be dealt with for EOR measurements (see Section 13.21 . 
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following two subsections we explore how sampling the 
statistical distribution can help separate the EOR sig- 
nal from the residual foreground and provide additional 
constraints on the process of reionization. 

3.2. Symmetries and Residual Foreground Subtraction 

The statistics of the model predictions for the EOR 
signal fall under the class of spatially homogeneous and 
isotropic random fields, since there is no preferred posi- 
tion or direction in space. The underlying isotropic sym- 
metry of the random field can be used to set strong con- 
straints on the expected signal properties. In the limit 
that the observed data are from the same epoch, the 
isotropy of the model implies that the Fourier transform 
representation of the neutral hydrogen emission will be 
spherically symmetric since the expectation is the same 
in all directions. The spherical symmetry of the EOR 
signal is an extraordinarily useful characteristic, and all 
the statistical characteristics of the signal will share this 
symmetry (including all moments of the distributions, 
not just the power spectrum distribution developed in the 
previous section, see Section [3. 3f) . The spherical symme- 
try can be used to differentiate the residual foreground 
signal from the EOR and to directly measure the statis- 
tical properties of the observed signal. 

Due to the spherical symmetry, all of the observations 
within a spherical shell of Fourier space are drawn from 
the same statistical ensemble. We can thus combine all 
of the measurements within a spherical shell in Fourier 
space to determine the underlying statistical distribu- 
tion. The spherical symmetry is only approximate be- 
cause the universe does evolve. However, the 4 MHz 
bandwidth of the baseline observations correspond to less 
than 4% change in redshift, over whi ch t he appr oxima- 
tion of no evolution holds quite well ifZaldarriag a et all 
|2004|) . Techniques which take into account the expected 
deviations from spherical symmetry may be employed in 
the future. 

Differentiating the very weak EOR signal from the 
foreground emission of the Milky Way and extragalac- 
tic sources is one of the primary difficulties with the 
EOR measurement. Planned EOR observations all in- 
volve selecting relatively cold and featureless portions of 
the Milky Way at high galactic latitude and carefully 
subtracting out the contributions of all known sources. 
However, the EOR signal is so weak that residual sub- 
traction errors and very fai nt foregroun d sources could 
still mask the desired signal l|Di Matteo et al. 2002). The 
question becomes how to separate the EOR signal from 
the residual foreground contamination. 

Fundamentally, the spherical symmetry of the EOR 
signal is based on the combination of narrow line emis- 
sion and redshift which allows one to form an isotropic 
image cube, where frequency maps to the line of sight 
distance (see Figure OJ. This mechanism is not present 
for most of the foreground contaminants, which either 
have smooth continuum spectra or are local line sources. 
As discussed in the following paragraphs, both of these 
foreground types have a separable- axial symmetry where 
the specific intensity Ip{u,v,ri) for the foreground has 
the form A{ri)B{\/u'^ + v^) (functions A and B are both 
complex). We can use the difference in symmetry be- 
tween the spherical EOR signal and the separable-axial 




Fig. 3. — Cartoon of the {I{u, v, r])) distribution of the EOR sig- 
nal. Notice the spherical symmetry of the EOR signal as compared 
to the strong separable-axial symmetry of the residual foreground 
sources as shown in Figures |4| and |F] This symmetry difference 
can be used as a handle for separating the EOR signal from the 
residual foreground contamination. 



residual foregrounds to distinguish the EOR signal from 
the foreground contaminants. 

The most problematic smooth spectrum foreground 
consists of faint emissi on from galaxies at cos mologi- 
cal d istances foi Matt eo et al.i r2002: Zaldarriaaa et all 
120041) . The spectra of these sources is smooth in fre- 
quency (no sharp lines) and the spectral power law and 
curvature characteristics are independent of look direc- 
tion. In the image cube, these sources map to all frequen- 
cies at a given spatial location, producing pencil lines of 
emission along the frequency axis (the contribution need 
not be uniform along the pencil beam). Spatial cluster- 
ing of these foreground sources can lead to features in the 
u and V Fourier coordinates, but not in the rj coordinate 
since the emission does not map to a single frequency. 
This creates the separable-axial symmetry in the Fourier 
space for the dominant emission from faint unsubtracted 
point sources (see Figure^. 

The strongest foreground signal is the diffuse syn- 
chrotron radiation from the Milky Way, though this 
foreground is probably one of the easiest to remove 
(Di Matteo et al. 2002). Because this emission is fairly 
smooth in both angle and frequency, the residual errors 
from this emission is concentrated at small wavenum- 
bers in the Fourier plane, as shown in Figure There 
may well be weak small scale structure due to clumps in 
the galactic emission and free-free absorption, and the 
correlation properties of this weak small scale structure 
is unknown. However, since the emission is broadband 
the resulting contamination has the same separable-axial 
symmetry in the Fourier representation as the smooth 
spectrum point sources. 

A more difficult local foreground is produced by radio 
recombination lines from our own galaxy. Line emission 
will also be included in the foreground model and sub- 
tracted, but significant errors in predicting and removing 
the recombination lines is expected. However, in the im- 
age cube the residual recombination lines will appear as 
planes at set frequencies (the magnitude and sign of a 
residual line may vary with position, but the frequency 
is set). This again leads to the separable-axial symmetry 
in the Fourier cube. 

Not all of the foregrounds have a separable-axial 
symmetry — any line source at cosmological distances will 
partially mimic the spherical symmetry of the EOR sig- 
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Fig. 4. — Cartoon of the {I{u, v, rj)) distribution of the resid- 
ual faint galaxy foreground. While there is significant power in 
the spatial directions u and v (Di Matteo et al. 2002), the spec- 
tral smoothness of the faint galaxy residual foreground leads to 
a strong separable-axial symmetry with the contribution concen- 
trated at small rj values. 




Fig. 5. — Cartoon of the {I{u,v,r])) distribution of the residual 
Milk Way synchrotron foreground. The contribution tends to low 
values because the residual Milky Way synchrotron foreground is 
smooth both spectrally and spatially, and has a separable-axial 
symmetry due to the spectral characteristics. 



nal. Fortunately, no strong lines at frequencies near or 
below the 21 cm line are expected in the foreground sig- 
nal. Sources such as radio recombination lines in high 
redshift galaxies will be a confusing source, but this emis- 
sion is expected to be very weak and will have a slightly 
elliptical symmetry due to differences in the frequency to 
redshift scaling. 

Additionally, the symmetry argument does not remove 
the issue of foregrounds. Large residual foreground sub- 
traction errors can swamp the EOR signal regardless of 
symmetry; the symmetry simply provides an additional 
handle for separating the foreground from the signal. 
The accuracy to which a given foreground needs to be 
subtracted depends on both its strength and the exact 
shape of the Fourier signal from that source. For exam- 
ple, unsubtracted discreet foreground galaxies are rel- 
atively bright, but because they are spectrally smooth 
their contribution is largely confined to small rj values. 
Radio recombination lines on the other hand are rela- 
tively weak, but due to the sharp spectral features may 
contain much more power at high r/ values, though this 
depends critically on the characteristics of the radio re- 
combination line subtraction process. 

The actual contribution of a residual foreground to the 
Fourier representation can be directly calculated given 
both a model for the foreground signal and the fore- 
ground subtraction process. Returning to the visibility 
cube V{u,v, Af), let us consider the errors made along 
the frequency dimension for fixed u and v. For smooth 



spectrum sources the two largest errors are going to be 
constant (complex) offsets at all frequencies due to resid- 
ual source subtraction errors, and residual spectral slopes 
due to errors determining the spectra of the foreground 
sources. This suggests a Taylor expansion of the residual 
errors along the frequency axis, with source subtraction 
errors providing the constant term, spectral index the 
first order term, and spectral curvature the higher order 
terms. The magnitude of the contamination is given by 
the probability of making a given error Pn(e), where e is 
the complex error and n is the order of the Taylor ex- 
pansion in the frequency dimension (in general Pn{t) can 
be a function of u and v). For the visibilities the error 
is given by sampling from the probability distribution 
P„(e) at each order 

oo 

V'{u, V, A/) = V{u, V, A/) + ^ P„(e) e A/", (23) 

n=0 

where V is the visibility from the EOR signal and V' in- 
cludes the residual foreground subtraction errors. Equa- 
tion can be Fourier transformed along the frequency 
axis to give the contamination observed in the Fourier 
representation 



Ai'iu, V, Tj) = Ai{u, v,r^) + Y, {Pn{e) e FT(Ar)} 

ji=0 



. (24) 

For many foreground sources and subtraction algorithms 
the first few terms of the Taylor expansion are expected 
to dominate. 

The spherical symmetry discussed here is largely in- 
dependent of the source subtraction method, and can 
be used in conjunction with other methods to further 
reduce the effect of residual foreground contamination. 
Because the proposed EOR observations are many or- 
ders of magnitude deeper than existing measurements, 
the brightness and spectral characteristics of many of the 
foreground sources are poorly constrained. This makes 
estimating the residual contamination left by particular 
subtraction algorithms very difficult. Measuring the fore- 
grounds and studying the residual subtraction errors is 
the subject of several ongoing research programs and pro- 
posed observing campaigns around the world. 

Our current understanding of the residual foregrounds 
can be summarized as follows. The smooth spectrum 
emission from the faint galaxy foreground can be re- 
moved using the spectral information, either by using the 
correlation of visibilities in widely separated frequency 
channels, or traditional spectral fitting techniques (Zal- 
darriaga et al. 2004; Briggs et al. private communication 
2004) . Both of these techniques will clean the faint point 
source foreground to a level less than the expected EOR 
signal, and the spherical symmetry can be used to re- 
duce it further. The smooth spectrum contribution from 
the Milky Way limits the sensitivity of the observations 
by being the dominant source of the thermal noise (see 
Section , however it is expected to be quite smooth 
and easily removed using traditional spectral fitting tech- 
niques. The residual foreground from radio recombina- 
tion lines is more problematic due to uncertainties in the 
emission characteristics. Radio recombination lines in 
our own galaxy will be measured and subtracted, but 
care must be taken to accurately determine the relative 
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strength of the different emission hnes to allow precision 
subtraction. Accurately measuring the Milky Way radio 
recombination lines is one of the main motivations be- 
hind several current observational campaigns. The radio 
recombination lines from faint galaxies must be removed 
to high accuracy because they can mimic the spherical 
symmetry of the EOR signal, however, they are quite 
weak. If radio recombination lines can be subtracted for 
sources above a few mJy and contribute 0.1% of the flux 
for sources weaker than this, their contribution to the 
total variance will be an order of magnitude fainter than 
the expected EOR signal. 

The strong symmetry differences between the expected 
signal and the residual foreground contaminants provide 
an additional handle for separating the EOR signal from 
most of the residual foregrounds. In practice, the sepa- 
ration is probably best performed by modeling the resid- 
ual effects of the foreground sources, and using a maxi- 
mum likelihood analysis to identify the EOR emission. In 
this type of analysis the difference in symmetry between 
the signal and the foregrounds provides one particularly 
powerful handle, which can be combined other statistical 
signatures to identify and measure the properties of the 
EOR. In future work we plan to extensively model the 
foreground contamination and removal process to fully 
explore the effect of foreground sources. 

3.3. Non-Gaussian Statistics 

Unlike the density fluctuations of the CMB, the EOR 
signal is expected to be non-Gaussian in nature. This 
is because the intensity of the EOR signal is a function 
of the spin temperature of the gas and the local ion- 
ization fraction, in addition to the density fluctuations. 
While the density fluctuations are largely Gaussian (at 
least in the linear gravity regime) , the spin temperature 
and ionization fraction may have a very complex distri- 
bution ( Madau. Mciksin & Rces 1997; Tozzi et al. 2000), 
and the overall statistical distribution is unlikely to have 
a Gaussian distribution. 

The power spectrum is still the Fourier transform of 
the autocorrelation function and as such is a uniquely 
important measure of the image statistics, and the de- 
velopment in Section is valid. But unlike with the 
Gaussian fluctuations of the CMB, the power spectrum 
of the EOR signa l does not fully desc ribe t he underly- 
ing statistics (see iCox &: Milled 11196 5') and ' Vanmarckd 
(p983) for background on the Weiner-Khinchine theorem 
and related issues for non- Gaussian random fields). In 
particular, the full shape of the statistical distributions 
within the spherical Fourier annuli — the power spectrum 
is just the second moment of these distributions — and 
the global statistical distribution in the image cube may 
provide important information on the physics of reion- 
ization, including the history of the spin temperature 
and the characteristics of the reionization front. At this 
point, there is little theoretical guidance on the shape 
that the statistical distributions may take, but it is im- 
portant to include these effects in future calculations as 
they may provide important information on the process 
of reionization. 

4. QUALITATIVE EFFECTS OF REIONIZATION ON THE 
POWER SPECTRUM 



A full prediction of the expected power spectrum and 
other statistical signatures of the EOR signal will re- 
quire detailed Monte Carlo simulations. However, there 
are a few general features of the EOR power spectrum 
which are expected. In this section we explore a couple of 
toy models to help build a conceptual understanding of 
the power spectrum characteristics and to demonstrate 
some of the information which can be extracted from the 
statistics of the EOR signal. 

The first signature is the power spectrum of matter 
density fluctuations. For a simple model with uniform 
fractional ionization and spin temperature, the EOR sig- 
nal is simply the primordial power spectrum (see Section 
Inland Tozzi et al. (2000)). Generally, the amplitude of 
this signal depends on the spin temperature, the frac- 
tional ionization, and the gravitational growth of density 
fluctuations. There are a number of effects which could 
also distort the observed shape of the density fluctua- 
tion power spectrum, but at large spatial scales (small 

k) local effects from non-uniform spin temperature and 
ionization bubbles should have little impact on the shape. 
Thus we expect the power spectrum of density fluctua- 
tions to be present at some amplitude in all models of 
the EOR (see Figure EJ. One enticing possibility raised 
by this signature is the opportunity to constrain the cos- 
mological equation of state if the effects of fractional ion- 
ization and spin temperature ca n be separated from the 
density fluctuation growth term iCnedin fc Shave'ill2004l: 
Loeb & Zaldarriaga 2004). 

The second signature is due to the Stromgren spheres 
which are expected to ionize the inter-galactic medium. 
These spheres are observed by the absence of neutral 
hydrogen emission on the ionized interior, and provide 
a strong power spectrum signal in most models de- 
veloped to date. Again, the the shape of small ion- 
ized regions may be strongly affected by local phenom- 
ena, but are expected to become somewhat spherical on 
large spatial scales as reionization progresses. The auto- 
correlation function of the Stromgren spheres will be 
strongly centralized, and the Fourier transform will have 
a sinc^ like distribution. Thus the process of reioniza- 
tion creates an increase in the power spectrum at small 
k with a shoulder at the characteristic length scale of 
the Stromgren spheres (see Figure Eland the simulations 
bv lZaldarriaea et al.l l)2004D ). As the spheres enlarge the 

shoulder will move to smaller k values and the amplitude 
of the power spectrum at small k will increase. Note also 
that the power spectrum will be distinctly different if the 
universe was ionized by a small number of quasars vs. a 
large number of luminous galaxies. If a small number of 
luminous sources dominate reionization, the ratio of the 
ionization fraction to the characteristic size of the ionized 
regions will be much lower than if a large number of less 
luminous sources dominate reionization. By measuring 
the relative amplitude of the the small k component of 
the power spectrum to the position of the shoulder, we 
can constrain how many sources participated in reioniza- 
tion of the universe (see Figure EI). 

These simple examples show the promise of studying 
the statis tical properties of the E OR signal. Very recent 
results bv Barkana fc Loebl ((2004) indicate that reioniza- 
tion was complicated and patchy, and the power on very 
large length scales is expected to increase significantly. 
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Fig. 6. — Two signatures which are expected to appear in the 
power spectrum of the EOR signal. The dashed hne is the power 
spectrum of cosmological density fluctuations (this example is at z 
= 10). While the shape of the distribution is robust, particularly at 
smaller baselines, the amplitude depends on the distribution of spin 
temperatures and the reionization history. The solid line is then 
a cartoon of the power spectrum produced by Stromgren spheres. 
The shoulder occurs at the characteristic size of the spheres and 
moves to the left as they become larger. At the same time the 
height of the low k portion will increase, and the ratio between 
these values can constrain the number of ionizing sources. 




Fig. 7. — The power spectrum of matter density fluctuations in 
observer coordinates u,v,t], with no ionization. The 21 cm emis- 
sion has been redshifted to 203 MHz, 142 MHz, and 109 MHz for z 
= 6, 9, 12 respectively. The horizontal axis has been scaled by the 
wavelength so that the visibility measurements made by antenna 
pairs remain constant in position for all three redshifts. In the cur- 
rent LOFAR design, the visibility measurements from the central 
core extend to distances of 2 km, though the distribution of these 
measurements has not been determined. 



Detailed Monte Carlo predictions of the expected statis- 
tical distributions are needed to allow even clearer differ- 
entiation between various models of structure formation 
and reionization. 

5. SIMPLEST CASE: REHEATED STRUCTURES BEFORE 
REIONIZATION 

In this section we use the analysis technique developed 
in Section |3| to determine the expected signal strength 
for a typical wide field radio observation. The model 
we will consider assumes a preheated spin temperature 
with no significant ionization and h as been extens ively 
discussed by Tozzi ct al. (2000), Ili ev et al.l l|2002|) and 
others. For the instrument we will use the baseline de- 
sign for the LOFAR observatory and assume a single 100 
hour wide field radio observation analyzed using only the 
power spectrum.^ This case is useful for comparing with 
pervious work and sets a baseline sensitivity for this anal- 
ysis technique. The exact sensitivity of any measurement 
depends intimately on the characteristics of the instru- 
ment. However, an estimate of the EOR sensitivity can 
be developed using a few approximations. 

For observations at 129 MHz (z=10) with a uniformly 
heated medium (Tg ^ Tcmb), ^hi is just the constant 
Cjy ~ 164 Jy^ str~^ times the matter density power 
spectrum (see Appendix ^ for a discussion of units). 
The power spectrum of the matter density fluctuations 
is shown in Figure indicating the strong signature of 
the density fluctuations which will be imprinted on the 
EOR signal. The function 77)^ from Equation 

^ T he full design of the LOFAR, observatory 
jhttpj/ /www. lofar.org I is currently in a state of flux, and it 
IS unclear whether the observatory being constructed in the 
Netherlands will be capable of wide field observations or have 
the full baseline sensitivity. However, this is a good baseline and 
future observatories are likely to be significantly more sensitive. 



1221 is very sharply peaked, so we can make the approx- 
imation that the power spectrum is a constant on the 
length scale of \W{u,v,ri)\'^ and pull the matter density 
power spectrum out of the integral. (Since in practice 
we are fully sampling the instrumental bandwidth, we 
are really performing a uniform average over 77 in this ar- 
gument and an equivalent answer can be obtained using 
a variation of Parseval's relation for 77.) If the response of 
the antennas can be described as constant to 20° with 4 
MHz bandwidth, the integral of \W{u, v, ry)^ is 3.8 x 10^ 
str Hz by Parseval's relation. The baseline LOFAR de- 
sign is centrally condensed and does not evenly sample 
the u,v,r] space, but the magnitude of the condensation 
has yet to be determined. If we assume that the sampling 
varies as approximately as where r 
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the average neutral hydrogen power spectrum is 0.74 str 
Hz. Combining all of these factors we obtain an average 

(\Ai{ui,vi,r]i)\^^ signal of 4.5 x 10^ Jy^ Hz^. 

The next step in calculating the sensitivity is to de- 
termine the thermal noise inherent in the measurements. 
The total thermal noise also depends on the distribution 
of antennas on the ground, but for a dense array may be 
approximated by the following argument. For a single 
dual polarization visibility measurement from a pair of 
high frequency antenna elements (each with 16 dipoles) 
the thermal noise is given by 



sys 



(25) 



where df is the width of the individual frequency chan- 
nels, and T is the duration of the observation. To calcu- 
late the noise on A/(u,w,7y), we need to Fourier trans- 
form AV{u,v, f). Since the noise in AI(u,v,ri) will be 
the same at all rj, we can calculate the noise at 77 = and 
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generalize to the other values. AI{u, v, 0) becomes 



A/(u,«,0)= / Viu,v,f)df. 



B/df 



Viu,vJ)df, 



(26a) 



(26b) 



where B is the fuU bandwidt h of th e measurement. The 
rms of the sum is ^/N = y/B/df times the rnis of V 
times df , and is given by 



A-^rnis — 



\/2kBTsysVB 



4A2 



(27) 



Using the values of T^ys = 695 K, A = 2.32 m, r = 100 
h, and B = 4 MHz gives A/,,„s = 2.1 x 10^ Jy Hz. 

The background noise {AI(u,v,r])) is Gaussian dis- 
tributed with rms a. The expectation of the absolute 
value (|AJ(m, v,ri)\) is then Rayleigh distributed and the 
expectation of the absolute value squared (| A/(it, v, ry)^) 
is exponentially distributed with an rms of 2a^. This 
gives an rms for a single |A/(w, u, ry)p measurement of 
8.8 X 10^" Jy^ Hz^. If we average over a large number of 

visibilities, the mean (| A/(m, u, 77)^) will become Gaus- 
sian distributed with width rms/ ^^N^is, where A'vis is the 
number of visibility measurements in the sum. 

For our example we have 1.2 billion visibility measure- 
ments from the wide field correlator, giving an rms of 
2.5 X 10^ Jy^ Hz^ for the thermal noise. Comparing to 
the expected signal strength of 4.5 x 10^ Jy^ Hz^ results 
in a signal to noise ratio of '--^ 15 for a single 100 h 4 MHz 
observation in the absence of significant foreground con- 
tamination, or equivalently a ~ Scr signal in 5 separate 
annular bins. This suggests that precision measurements 
of the power spectrum can be achieved using wide field 
radio observations and reasonable integrations. 

6. CONCLUSION 

In this paper we discussed the power of wide field radio 
observations and introduced an analysis technique which 
we expect will efficiently extract the statistical properties 
of the EOR signal and separate them from the residual 
foreground contamination. The next step in developing 
the experimental procedure is to use simulations to ex- 
plore a number of subtle experimental effects and their 
impact on EOR observations. Using end-to-end simula- 
tions which include realistic foreground models and in- 
strumental effects, we intend to explore four major top- 
ics: 

• Foreground contamination and our ability to mea- 
sure the EOR power spectrum using various fore- 
ground subtraction methods and symmetry tech- 
niques. 

• Ionospheric effects and calibration errors, and the 
impact on deep EOR integrations. 

• Array optimization to enhance the science return 
of EOR measurements in future observatories. 

• Statistical properties of theroetical models of reion- 
ization, and the ability of wide field radio observa- 
tions to differentiate between EOR models. 



It is clear that wide field radio observations can make an 
important contribution to the study of the early universe, 
and the above steps will help formulate a robust strategy 
for measuring and interpreting the EOR signal. 
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A. DERIVATION OF BRIGHTNESS TEMPERATURE 

The optical depth is related to the absorption co- 
efficicnt by (see, for example, iRvbicki fc LightmanI 
^979,) Equation 1.26) 



Tuis) = / a^{s')ds'. 

J So 



(A-1) 



Under conditions of local thermodynamic equilibrium ap- 
plicable for a thermalized neutral hydrogen gas the ab- 
sorption coefficient is given by 



= (1 



-hu/kTs 



(A-2) 



from lRvbicki fc LightmanI l)1979|) Equation 1.80, where v 
is the frequency, ni the number density of lower quantum 
states, B12 the Einstein absorption coefficient, T, is the 
temperature of the spin states of the neutral hydrogen 
and (f){v) is the line profile. The relation between the 
absorption coefficient B12 and the transition probability 
for spontaneous emission A21 with the number of states 
in the upper and lower energy levels given by 32 and gi 
respectively is 



2Ju^ 
A21 — — ^ i>i2, 



-hv/kTs 



(A-3) 



(A-4) 



which we can then use to obtain 

92 . (. 
= 7, — 5-— "^1^21 1 
87rz^^ g\ V 

For the 21 cm line of hydrogen, the factor hv jk is 0.068 
K. This is a much lower temperature than the CMB ra- 
diation, and for almost all cases the spin temperature 
will be warmer than this even when the 21 cm line is seen 
in absorption. Making the approximation hv /k <^Ts 'we 
find 

92 A 

Using the same approximation we can also show that the 
hydrogen will be equally distributed between all of the 
Si quantum states. For the 21 cm transition 52 = 3 and 
gi — I. This then implies that the ni = nHi/4, where rifji 
is the number density of neutral hydrogen at the given 
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redshift and implicity includes the density fluctuations 
(the units have not been defined yet). Inserting these 
values we obtain 

Now to calculate the opti cal de pth we need to perform 
the integration in Equation I A- II Because the 21 cm line 
is extremely narrow all of the optical depth is constrained 
to a specific redshift z. Effectively (f){v) becomes a 5- 
function over the line of sight variable s' . We can thus 
rewrite the equation using 4>[v) = 5{v2i ^ (1 + z')v) = 
^^Xl ~ '^'•^ ^^'^ perform the integral to obtain 

"ic? ^ hv-ix (1 + z') ^ , ..ds , . 



3c2 



327ri^2i 



nmA 



21 



kTs V2i 
hv2i (1 + z) ds 
kTs 1^21 dz' 



(A- 



where we have evaluated the frequency terms from the 
transition strength at the 21 cm line and left the ds/dz 
conversion from length sc ale to z until we explicitly de- 
cide on units in Equation lA- 1 21 below. 

We can now use the optical depth to obtain the spe- 
cific intensity at the s ource of the 21 cm emission. From 
IRvbicki fc LightmanI ljl979(l Equation 1.30 the specific 
intensity after the emission region (ly) is given in terms 
of the background brightness /cmb, which is dominated 
by the cosmic microwave background, and the emission 
brightness of the neutral hydrogen /hi via 



-^CMBe 



/hi 1 



(A-9) 



The difference between the background and the neutral 
hydrogen emission is then given by 



A/„ = (/hi - /cmb) (1 



(/] 



HI 



/cmb) Ti/, 

(A-10) 

where the last relation follows from the small optical 
depths associated with the 21 cm line. 

The final step in the process is to propagate the radia- 
tion from the emission region to the Earth. The number 
density nm is related to the mean hydrogen density (nn) 
times the neutral hydrogen density fluctuations 



"-HI 



Pm 

(ph)' 



(A-11) 



We now have to decide on the units we will use for nni, 
ds/dz and the relativistic propagation relationship. Here 
we choose the units of the emission region to be proper 
Mpc, though all unit choices will give the same answer. 
Since the average neutral hydrogen density is locked to 
the Hubble flow, the mean density evolves as (1 -I- zY . 
The mean hydrogen density can then be defined in terms 
of the present value of the cosmological baryon density 
f^B to obtain 

(nn) = P^{l~Y){l+zf = {n^h^) (l-r)|^/^(l+z; 

where Y is the helium fraction, G is the gravitational 
constant, ps is the physical baryon density, H^/h is the 
Hubble constant, and mn is the mass of a hydrogen atom. 



The proper length per redshift interval ds/dz is given by 



ds 
dz 



^ (A-13) 

yJVLM{l + zf + f^k(l + 2)2 4- VLa. In 



H^{\+z)FAzy 



where E{z) 

proper coordi nates the brigh tness conservation law is 
given by ffrom lPeacocS l|1999D Equation 3.88) 



/.(O) 



(1 + Z)3 



(A-14) 



Putting together Equations IA-81 and lA-lOl through IA-141 
and simplifying one then obtains 



3c^ 



327ri/^i 



KHo/hf 



Btt G mn 



A21 hv2i (/hi — /cmb) Pm 



Ho E{z) k 



(A-15) 



The specific intensities can be replaced by brightness 
temperatures to obtain 



AT, = 



3c^ 



Ml hv2i {Ts - Tc 



MB I 



HHo/hT 

Svr G mn 
2 Phi 



Ho E{z) k 



(1 + ^) 



(ph)' 



(A-16) 



where T, is the spin temperature and the (1 -I- z)^ factor 
relates the temperature on the left hand side which is 
measured at the antenna and the temperatures on the 
right hand side which are in the frame of the neutral 
hydrogen. 

Substituting standard values for the constants [Vl-Qh? 
= 0.02, Y = 0.24) these simplify to 



A/, = (2.9 mK) h 



AT, = (2.9 niK) h 



-1 ( /hi — /cmb) Phi 
E{z)Ts (ph)' 



(A-17) 



_i (1 -I- zf (Ts - Tcmb) Phi 



E{z) 



(ph) 



(A-18) 



B. UNITS 



The measured surface brightness in sky coordinates 
AI{6x, 9y, f) is in Jy str~^. The window function in sky 
coordinates W{9x,9y, f) is unitless. The Fourier trans- 
formed surface brightness AI{u,v,ri) is in Jy Hz, and 
the window function W{u,v,ri) is in str Hz. For equa- 
tions such as number 1151 the convolution over du dv drj 
removes the extra factor of str Hz (see Footnote QJ. The 
measured visibility V{u,v,f) is measured in Jy. 

The quantities Cjy and Ck are measured in Jy str~^ 
and K respectively. The neutral hydrogen density 
^"(ph)*"^ unitless in 9x,6y,Af coordinates, and has 
3 units of str Hz in Fourier coordinates. The delta func- 
' tion S{n — n'") has units of str~^ Hz"^, leading to the 
power spectrum PHi(fi) having units of Jy^ Hz str^^ or 
str Hz. If the spin temperature is a constant, C can be 
removed and the power spectrum of the neutral hydrogen 
density fluctuations is str Hz as one would expect. 
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